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Abstract 

A study of secondary Drell-Yan production in nuclear collisions is presented 
for SPS energies. In addition to the lepton pairs produced in the initial 
collisions of the projectile and target nucleons, we consider the potentially 
high dilepton yield from hard valence antiquarks in produced mesons and 
antibaryons. We calculate the secondary Drell-Yan contributions taking the 
collision spectrum of hadrons from the microscopic model URQMD. The con¬ 
tributions from meson-baryon interactions, small in hadron-nucleus interac¬ 
tions, are found to be substantial in nucleus-nucleus collisions at low dilepton 
masses. Preresonance collisions of partons may further increase the yields. 
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I. INTRODUCTION 


A unique experimental signature from the early stage of violent heavy ion collisions is 
provided by dilepton radiation. The (thermal) radiation of leptons is highest when the en¬ 
ergy density has reached a maximum. Since leptons are insensitive to strong interactions, 
the signal is not distorted during the further evolution of the system. Recent reports on in¬ 
termediate and high mass muon pairs produced in heavy ion collisions have attracted 
much attention. The suppression of the signal, as extracted from NA38/NA50 measure¬ 
ments, might indicate that deconhnement has been reached in central Pb-|-Pb reactions [^. 
In the mass region below the p-meson, on the other hand, the measured e'''e“ pair spectrum 
is enhanced relative to the expected ‘cocktail’ of hadronic decays @]. A good understanding 
of the continuum dilepton spectrum, in the relevant invariant mass window, is a prerequisite 
for the correct interpretation of the data. 

Rate estimates of Drell-Yan pair production have been performed systematically for 
different systems, energies, and invariant masses 0. They are essentially based on ex¬ 
trapolations of p-|-A reactions where a linear scaling with A is observed for the Drell-Yan 
production cross section [^: 

= o'o A . 

This linear scaling can be understood in the Glauber picture of the hadron-nucleus cross sec¬ 
tion, constructed at high energies using the AGK cutting rules . (For a simple dynamical 
picture of this effect, see |^.) 

A number of different theoretical frameworks have been used to explore the degree 

to which the hot matter produced in ultrarelativistic heavy ion collisions contributes to the 
total dilepton radiation. Dilepton emission from thermal sources such as a quark-gluon 
plasma or hadron gas, whether in or out of equilibrium, cannot be neglected relative to 
nucleon-nucleon contributions to the Drell-Yan process for masses below the J/V' peak — 
at least at higher energies than those presently available. Another source of dileptons which 
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has not been fully explored is Drell-Yan production by interactions between produced, or 
secondary, hadrons. Unless these secondary dilepton sources are understood, the thermal 
contributions cannot be reliably extracted. In this paper, we study possible rehnements 
to Drell-Yan production in nucleus-nucleus collisions from secondary hadron interactions at 
current hxed-target energies. It is shown that these secondary collisions can serve as an 
important source of m ~ 2 GeV dileptons due to the availability of valence antiquarks in 
mesons and antibaryons. 


II. THE MODEL 


In order to investigate secondary dilepton production at SPS energies^] a microscopic 
hadronic transport code, URQMD [^, is employed to obtain a realistic collision spectrum 
of secondary hadrons. The differential Drell-Yan cross section is computed at leading order 
(LO) using the standard equation ||T7|: 


^ ^ -(AB ^ IIX) = el \q^{xA, m^) + q^{xA, m‘^)q^{xB, m^) 


Qrn?s q ^ 


( 1 ) 


dw?dy 

where g(x, m?) and q{x, m?) denote the quark and antiquark densities, ^/s is the center- 
of-mass energy of the colliding hadrons, m is the invariant mass of the lepton pair, xa = 
\/r and xb = a/t with r = m? js, and y is the dilepton rapidity in the cms frame. 
The convolution of the parton distributions is weighted by the square of the charge of the 
annihilating quarks, Cq. Since the calculation is performed at leading order only, the result is 
multiplied by a “K-factor” [|T^. Although the ratio between the calculated next-to-leading 
order (NLO) and LO Drell-Yan cross sections is a function of m? and y, the dependence is 


^At RHIC and LHC energies the heavy flavour (c and b) combinatorial decays are expected to 


dominate the dilepton continuum |15|. Moreover, the nuclear parton distribution functions at such 
small X values are still quite uncertain. Therefore, we will not report any high energy calculations 
in this paper. 
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small, ~ 25% for m > 1.5 GeV and 0 < ?/ < 1 at a/s = 19.4 GeV. Thus, for simplicity, we 
have used a constant K factor. This is reasonable considering the accuracy of the present 
experimental data which does not reveal any strong mass dependence of the experimental 
K factor [| (the ratio between the measured data and the LO calculation). 

According to the URQMD hadronic model, most of the interacting particles in the early 
stage of the reaction are baryon and meson resonances. In this study, if two hadrons have the 
same quark content, their parton distributions are assumed to be equivalent, e.g. the proton 
parton distributions are used for the A^*(1440) and the pion parton distributions are used for 
the p. We have employed the GRV LO parametrizations of the parton distribution functions 
of nucleons and pions [^, since they provide a consistent set for both species (same 
Aqcd and similar initial scale /iq) . The lepton pair production cross section is calculated for 
each hadron-hadron collision and weighted by the inverse of the total hadron-hadron cross 
section. The distributions from these elementary hh-collisions are then summed. 

It is clear that in pion-nucleon, pion-pion and antinucleon-nucleon collisions, valence 
quark-antiquark annihilation can play a signihcant role in the Drell-Yan process. Figure 
shows the parton distribution functions of protons and pions from the GRV LO 
parametrization, for two different values of m?. For x ~ 0.05 the probability to hnd an 
antiquark is much higher in a pion than in a nucleon. The pion-nucleon cross sections 
are consequently higher than the nucleon-nucleon cross sections, especially when m/^/s ~ 
0.1. (See 1^.) On the other hand, in heavy ion collisions the typical meson-baryon and 
antibaryon-baryon interactions are certainly less energetic than the primary nucleon-nucleon 
collisions. 

While the usual Drell-Yan calculations of nucleus-nucleus collisions only consider in¬ 
teractions of the initial nucleons, in the present scenario four different processes may con- 


^This situation might change in the near future when the very high statistics Drell-Yan differential 
cross sections of E866 |[^ become available. 
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tribute: baryon-baryon, meson-baryon (and meson-antibaryon), meson-meson ^andbaryon- 
antibaryon interactions. Thus, it is interesting to estimate the relative importance of these 
contributions to the total Drell-Yan cross section. 

The typical scaling of the Drell-Yan cross section with target mass observed in hadron- 
nucleus collisions, accounted for here, assumes an equal probability for Drell-Yan pair pro¬ 
duction in each nucleon-nucleon interaction. In addition, in a separate microscopic simu¬ 
lation the meson-baryon, meson-meson and baryon-antibaryon collisions, involving newly 
produced hadrons, are calculated within the URQMD transport model. In order to conserve 
energy and momentum, the energy loss of baryons which have already interacted is taken 
into account in these secondary interactions. 

For simplicity, the produced hadrons follow straight-line trajectories across the reaction 
zone without energy loss (0°-scattering, only elastic interactions). The center of mass energy 
of these collisions is taken directly from the fully dynamical evolution of all hadrons within 
the URQMD simulation. In the usual cascade prescription, leading hadrons, containing 
valence quarks from the incident nucleon, can interact within their formation time with 
a reduced but finite cross section. Since the constituent antiquark content is much more 
relevant for the Drell-Yan process, we only consider collisions of formed hadrons. The 
formation time of produced hadrons, on average Tp ^ 1 fm/c, is determined by the string 
fragmentation prescription used in URQMD (defined as the hrst crossing point of the 
quark and antiquark or quark and diquark which form the hadron). The exact formation 
time depends on the mass and momentum of the produced hadron. By imposing a minimum 
hadron formation time, we effectively neglect any possible annihilation of produced partons 


^The Drell-Yan process in meson-meson interactions should not be confused with processes such 
as 7r+7r“ ^ l~^l~ where the whole meson is involved in the annihilation rather than a single parton. 
In any case, dileptons from pion annihilations are not expected to make an important contribution 
to the intermediate mass region we are concerned with here. 
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which have not yet hadronized (“preresonance” scattering). In Sec. ^ we return to this 
problem. 


III. THE EFFECT OF A FLAVOUR ASYMMETRIC SEA 

Important consequences for the sea quark distributions have recently been discussed: 
modihcations at small x and the violation of the Gottfried sum rule |^. Two effects are 
primarily important at low x. The HERA data has shown a steep rise in the proton structure 
function F 2 for x ~ 0.01 Additionally, the sea quark and gluon distributions can be 
signihcantly modihed in nuclear targets The observed low x shadowing and the low x 
rise in the parton sea distributions are not relevant for SPS energies where 0.077 < x < 0.52 
for 1.5 < m < 10 GeV at midrapidity. At higher energies, such as those at RHIG and 
LHG, these effects cannot be neglected. However, as we now describe, the violation of the 
Gottfried sum rule has important consequences at the SPS energies. 

The Gottfried sum rule dehnes the difference between the proton and neutron structure 
functions [E2l 


Fj") 


losR = E(0,1) = /" —(Ff 

Jo X 

IF 2 F - 1 - 

= - dx{uv — dv) H— / dx{u — d) = - if u = d. 
3 Jo 3 Jo 3 


( 2 ) 


Since the neutron distributions are not directly measured, is obtained assuming that 
M^(x,m^) = d^{x,m‘^) and uP{x,rn?) = d"'(x,m^). All global structure function analy¬ 
ses prior to 1992 assumed that the light-quark sea distributions are flavour independent, 
u{x,m?) = d{x,m‘^). On the contrary, based on their measurements, NMG found 

that 

/* 0.8 

E(0.004, 0.8) = / —(F| - F^) = 0.227 ± 0.007 (stat.) ± 0.014 (syst.) (3) 

2 o .004 X 

at m? = 4 GeV^. A straightforward comparison implies that d > u. These results have been 


conhrmed by p-|-p and p-|-d asymmetry studies of Drell-Yan production by NA51 and 
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by preliminary results from E866 . The lack of a Regge f — A 2 exchange degeneracy also 

leads one to expect that 


A = d — u ^ X 


(4) 


at small x, where the Regge intercept is aji ~ 0.5. 

With this in mind, we compare two versions of the GRV proton parametrizations, the 

and the 1994 version that explicitly takes into 


flavour-symmetric set from 1992 


account the NA51 results by choosing 

M = ^ {u + d) - A 


d — — (u -|- d) -|- A 


(5) 


Another factor contributing to the difference between the two sets is the starting distribution 


used to obtain the two GRV valence distributions: the flavour symmetric KMRS B- [^] for 
GRV 92 and the flavour asymmetric MRS A [Q for GRV 94. These starting distributions 
differ at low x and in the fitted value of Aqcd as well as in their treatment of the u and d 
distributions. To study the difference between these two GRV parametrizations, we evaluate 
the lepton pair production cross sections for p-|-p, p-|-n and n-|-n interactions. (Note that at 
y = 0, p-|-n = n-|-p.) The ratios of cross sections from the GRV 94 relative to the GRV 92 
parametrizations are shown in Fig. ^as functions of m at ?/ = 0 and ^/s = 15 GeV. The effect 
of the flavour asymmetric sea, most pronounced for p-|-n and n-|-n interactions, is clearly 
important. Inclusion of the flavour asymmetric sea decreases the isospin dependence of the 
Drell-Yan cross section. The isospin correction needed for Pb-|-Pb interactions is 1.3 with 
the GRV 92 LO set 0] and 1.03 with the GRV 94 LO set. The smaller isospin correction is 


also found when MRS G distributions 31 are used 


Since neither corresponding d—u data nor parton distributions with a flavour asymmetric 
sea exist for pions, we use the flavour symmetric GRV distributions for the pion 


m our 


meson interactions. One might suppose the asymmetry to be less pronounced for the pion 
since the n^lud) and 7i~{ud) each contain a u and a d quark or antiquark. 
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IV. PROTON-NUCLEUS AND PION-NUCLEUS REACTIONS 


Figure ^ shows the mass spectrum of dimuons in p+Pt collisions at piab = 400 GeV/c for 
three different longitudinal momentum windows. Our calculation with the GRV 94 parton 


distributions provides a reasonable description of the NAS data [^| with K = 1.5. In addi¬ 
tion we have compared the model calculation with p+Gu collisions at piab = 800 GeV/c 
and found equally good agreement. In fact, the Drell-Yan mass and xp double differential 
cross sections measured by NAS [^, E605 and E772 experiments, are all in good 
agreement with each other, as shown in Ref. P5[] . 

The proton-nucleus calculations reveal that the contribution from secondary collisions is 
less than 10% of the dilepton yield for masses m > 1.5 GeV. The relative importance of 
secondary contributions decreases with increasing xp- At xp = 0.S25 and xp = 0.625, the 
original Glauber type calculation is increased by less than 1 %, even at the lowest masses. 

The dimuon spectrum produced in tt” + IF reactions, at piab = 125 GeV/c, as measured 
by the E537 experiment [Q, is presented in Fig. ^ for a central and a forward xp bin. The 
results of our leading order calculations, scaled by a K factor of 2, are shown. A different 
K factor relative to the value used in the proton-induced reactions is needed since we use 
the GRV 92 pion set with the GRV 94 proton set. Although the measurements have large 
uncertainties, the agreement between our calculations and the data supports the validity of 
our treatment of secondary meson-baryon collisions. 


V. NUCLEUS-NUCLEUS COLLISIONS 

Figure shows the relative importance of the secondary dilepton production, compared 
to nucleon-nucleon collisions alone, for different projectile/target combinations. In the mass 
region around 2 GeV the additional contributions are not negligible in nucleus-nucleus col¬ 
lisions while the increase is less than 5 % in the p+W case. For masses above 3 GeV, the 
additional production beyond nucleon-nucleon collisions from our hadronic calculation is less 






than 10 % in both S+U and Pb+Pb interactions. Evidently the high mass dilepton yield 
is not signihcantly enhanced by the relatively low energy secondary scatterings. At masses 
below 2 GeV, on the contrary, it appears that the contribntion from collisions of newly pro- 
dnced particles cannot be neglected. At masses aronnd 1.5 GeV, an enhancement of 25% is 
expected in S+U interactions and 45 % in Pb+Pb interactions. Glearly, the importance of 
secondary contribntions to the Drell-Yan cross section increases with the system size. Note, 
however, that this enhancement is insnfficient to acconnt for the discrepancy reported in 
Ref. m between the data and the “expected sonrces”. 

Althongh the pertnrbative QGD calcnlation of dilepton prodnction at 2 GeV and below 
can only be used as a qualitative guideline, it is clear that the low mass dilepton continuum 
in nuclear collisions can be enhanced relative to a linear superposition of nucleon-nucleon 
interactions. At least part of the observed enhancement of muon pairs at intermediate and 
low masses |l[] might be caused by this previously neglected hadronic source. The perturba¬ 
tive calculation is dehnitely unreliable for masses below 1 GeV where a substantial increase 
in electron pair production in nucleus-nucleus collisions has been observed |^. To emphasize 
the theoretical uncertainty of the perturbative calculation for masses below 2 GeV, in Fig. 
and in the following hgures, we use thinner lines in that region. 

The relative contributions to the calculated dilepton mass spectra at ycms = 0.5 in central 
S+U and Pb+Pb collisions are presented in Fig. p. The decomposition of the secondary 
dilepton yield into each channel (meson-baryon, meson-meson, baryon-antibaryon) shows 
that by far the most important secondary contribution is from meson-baryon interactions. 
For masses above ~ 3 GeV they are virtually the only source of secondary dileptons. At 
midrapidity and for masses ~ 1 GeV, they account for at least 90 % of the secondary yield 
while meson-meson collisions contribute ~ 10 % and baryon-antibaryon collisions contribute 
less than 1 %. 
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VI. PRIMORDIAL SOURCES OF LEPTON PAIRS AND FORMATION TIME 


The standard Drell-Yan process corresponds to the interaction of fnlly formed hadrons. 
However, it was shown 0 that, dnring the early stages of the system evolntion, partons 
can scatter and annihilate before they have come on mass-shell Thns, we postnlate the 
existence of another dilepton sonrce: the annihilation of qnarks and antiqnarks which are 
not yet bonnd in a color-singlet hadron. 

To estimate the importance of these “primordial” or “preresonance” qq annihilations, 
we have calcnlated the contribntion of snch processes assnming that the asymptotic parton 
distribntion fnnctions are also valid for the primordial states. This assnmption is distinct 
from thermal dilepton prodnction where the quark and antiquark are in a thermal environ¬ 
ment, in or out of equilibrium, with temperature-dependent parton densities. We thus relax 
the restriction that the partons can only interact after they have hadronized. This is done 
very simply in our microscopic cascade calculation by decreasing the formation time of the 
produced hadrons, tf, from the “default” value of around 1 fm/c. 

The importance of this primordial (preresonance) contribution to the dilepton mass 
spectra is shown in Fig. 0. We see that when we set tf = 0, the secondary sources of 
dileptons become much more important than previously indicated. Indeed, the secondary 
dilepton yield increases by a factor of ~ 5 at all masses compared to the calculations with 
the default tf ^ 1 fm/c. These primordial interactions also affect p-|-A reactions, enhancing 
the secondary contribution. For example, the dilepton yield at m = 4.5 GeV and xf = 0.025 
shown in Fig. increases by 25 % when tf = 0. Nevertheless, good agreement with the NAS 
data is retained. 

In Fig. ^ we have also included a calculation with an intermediate formation time for 
the produced hadrons, tf = 0.5 fm/c, in order to study the sensitivity of the calculations 


^There are some similarities between the off-shell interactions of the primordial partons and the 


soft parton annihilation model proposed to explain the “anomalous” low mass lepton pairs |37|. 
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to this (model dependent) parameter. As we might expect, the reduced average formation 
time leads to an increased dilepton yield due to the larger number of possible hadron- 
hadron interactions. With this value of the formation time, the enhancement in the range 
1.5 < m < 2 GeV shows quite good agreement with the intermediate mass data |^. On the 
other hand, an inhnite formation time, Tp = oo, corresponds to the exclusion of all secondary 
interactions, in agreement with the Glauber-type Drell-Yan expectation. A recent work 
has already pointed out the strong influence of the hadron formation time on dilepton 
production via 7r’''7r“ annihilation. 


VII. CONCLUSION 

Our study indicates that the lepton pair continuum cross section for masses up to 3 GeV 
should not be interpreted solely in terms of Drell-Yan production in nucleon-nucleon inter¬ 
actions, even at SPS energies. The yield from secondary contributions depends strongly 
on the rapidity region and, of course, on the dilepton mass. Preresonance interactions are 
estimated to enhance the secondary yield by up to a factor of 5. 

Because meson-baryon interactions dominate secondary lepton pair production, the rela¬ 
tive importance of secondary dilepton sources evidently increases linearly with the number of 
produced pions at these energies. It may thus be possible to distinguish the secondary yield 
from thermal dilepton sources which are expected to scale quadratically with the pion mul¬ 
tiplicity. The multiplicity dependence also indicates that the dilepton yield from scattering 
of secondaries is likely to be considerably higher at RHIG and LHG energies. 
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FIGURES 



X 

FIG. 1. The GRV LO parton distribution functions of protons and pions for valence and 
sea quarks. We show the values for (u + d) quarks in the proton (dotted) and in the pion (full lines), 
and the values for u quarks in the proton (dashed lines). The thick lines indicate the distributions 
for m? ~ 0.3 GeV^, the thin lines those for = 9GeV^. Note the hard valence antiquark 
distribution in the pion. The valence antiquark distributions in the antibaryon correspond to the 
valence quark distributions in the baryons. 
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FIG. 2. Drell-Yan cross section calculated with the recent GRV 94 LO parametrizations |2C] 
where u / d, relative to the calculation done with the GRV 92 LO set which assumes a 
symmetric sea. Notice that the K factor to be applied to these calculations is 2 for GRV LO 92 


and 1.5 to GRV LO 94. 
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FIG. 3. Dilepton mass distributions from p+Pt collisions at piab = 400 GeV/c calculated for 
three xp values. The contribution of baryon-baryon collisions is shown (solid line) separately from 
the total spectra (dashed line) which includes secondaries. The data points are from experiment 
NAS 1^. We use a K factor of 1.5 for NN collisions (calculated with the GRV 94 nucleon PDFs) 
and of 2 for pion induced reactions (calculated with the GRV 92 pion PDFs). We use thinner lines 
for m < 2 GeV to emphasize the relatively large theoretical uncertainty in the calculation. 


18 







FIG. 4. Mass spectrum of dileptons for two xp bins in 7r“ + VF collisions at piab = 125 GeV/c. 
We show our leading order calculations with data from experiment FNAL-537 We use a K 
factor of 2 for these pion induced reactions. 
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a) 



FIG. 6. Differential Drell-Yan cross section, 
central collisions of S+U at piab = 200 A GeV 
The different Drell-Yan contributions from the 


b) 



per unit rapidity, at i/cms = 0.5, for the 10 % most 
(left) and of Pb+Pb at piab = 160 A GeV (right), 
hadronic simulation are shown. 
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FIG. 7. Ratio between the total dilepton mass spectrum and the nucleon-nucleon scattering 
alone, at ycms = 0.5, for S-l-U colisions at pi^h = 200 A GeV (left) and for Pb+Pb collisions 
at piab = 160 A GeV (right). The calculations correspond to different formation times for the 
secondary hadrons. 
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